Homeostatic competition drives tumor growth and metastasis nucleation 
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We propose a mechanism for tumor growth emphasizing the role of homeostatic regulation and 
tissue stability. We show that competition between surface and bulk effects leads to the existence of 
a critical size that must be overcome by metastases to reach macroscopic sizes. This property can 
qualitatively explain the observed size distributions of metastases, while size-independent growth 
rates cannot account for clinical and experimental data. In addition, it potentially explains the 
observed preferential growth of metastases on tissue surfaces and membranes such as the pleural 
and peritoneal layers, suggests a mechanism underlying the seed and soil hypothesis introduced by 
Stephen Paget in 1889 and yields realistic values for metastatic inefhciency. We propose a number 
of key experiments to test these concepts. The homeostatic pressure as introduced in this work 
could constitute a quantitative, experimentally accessible measure for the metastatic potential of 
early malignant growths. 



Introduction 

The progression of cancer is a multi-step process. Over 
80% of malignant tumors are carcinomas that originate 
in epithelial tissues from where they invade through the 
basal membrane into the connective tissue. At some 
point, subpopulations of cells may detach from the pri- 
mary tumor and spread via the bloodstream and the lym- 
phatic system. Some of them give rise to metastases 
in distant organs. Metastases account for the major- 
ity of patients' deaths due to cancer, and thus under- 
standing the metastatic process is of critical importance. 
The metastatic cascade is a very inefficient process, as 
only one in about a thousand cells that leave the pri- 
mary tumor goes on to form a macroscopic secondary 
tumor. This property is referred to as "metastatic inef- 
ficiency" (Chambers et al, 2002; Sahai, 2007). Recent 
experimental results have shown however that cell ex- 
travasation is highly efficient, namely that over 80% of 
the metastatic cells that are present in the bloodstream 
manage to enter a distant organ (Luzzi et al., 1998; 
Cameron et al., 2000; Zijlstra et al., 2002). Thus, the 
main contribution to metastatic inefficiency arises from 
the failure of cancerous cells to grow inside invaded or- 
gans. Metastatic tumors also show preferential growth 
in different organs with a distribution that cannot be 
explained by blood flow patterns alone. Hence, the ef- 
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ficiency of the metastatic process depends on specific 
interactions between the invading cancer cells and the 
local organ tissues (Fidler et al., 2003). This concept, 
referred to as the "seed and soil hypothesis" , was in- 
troduced by Stephen Paget as early as 1889 (Weinberg, 
2007) : "the seed" — the metastatic cell — needs to be com- 
patible with "the soil" — the host tissue — for successful 
growth to occur (Fidler, 2003; Couzin, 2003). Despite 
this early observation, the nature of the interactions con- 
trolling both the efficiency of the metastatic process and 
its tissue specificity remains a poorly understood aspect 
of cancer progression even today. 

In this work, we introduce the notion of homeostatic 
pressure and propose that it is an important property 
for describing the competition between different tissues 
grown in a finite volume. The concept of homeostatic 
pressure is best defined from the following experiment: 
consider a chamber in which cells can be cultured with 
a setup that enables successful proliferation, allowing in 
particular for water, oxygen, nutrients and growth factors 
to diffuse through the compartment walls, keeping the 
cells' biochemical environment constant. The compart- 
ment is closed on one side by a piston connected to a rigid 
wall with a spring (Fig. [iJV). As the growing tissue fills 
the available space and gradually compresses the spring, 
the pressure rises until a steady state is reached in which 
division balances apoptosis and the piston stops mov- 
ing. The spring position is stable, since further growth 
increases the pressure above this value and favors apopto- 
sis, whereas recession of the piston decreases the pressure 
and favors division. This steady state is characterized by 
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FIG. 1: (A) Schematic representation of a measurement ap- 
paratus for the homeostatic pressure. As the tissue prolifer- 
ates, the piston compresses the spring and the pressure ex- 
erted on the tissue increases. Once a steady state is reached, 
cell division and apoptosis balance. The cell density and the 
pressure exerted on the spring at this point define the home- 
ostatic mechanical state of the tissue in a given biochemi- 
cal environment. (B) Schematic representation of a tissue- 
competition experiment. The two tissues are in mechanical 
contact through a freely-moving, impermeable piston. The 
tissue with the lower homeostatic pressure is compressed to 
a cell density above its homeostatic point and initiates apop- 
tosis. The other tissue proliferates and expands until the op- 
posing tissue has disappeared. 



a well-defined pressure exerted on the spring and a well- 
defined density of cells, which we refer to as the homeo- 
static pressure and density of the tissue in this particu- 
lar biochemical environment. Note that the homeostatic 
pressure is different from the hydrostatic pressure since 
the lateral walls of the chamber allow for fluid transport. 
Instead, it resembles more an osmotic pressure but orig- 
inates from the forces driving tissue expansion. 

We now show that the ability of one tissue to replace 
another one in a competition for space depends on the 
relative values of their homeostatic pressures. Let us con- 
sider a similar chamber, but in which the piston separates 
two tissue compartments H and T, establishing mechan- 
ical contact (see Fig. [l]3). Suppose that tissue H has a 
homeostatic pressure PH,h smaller than that of tissue T, 
PT.h- As cells divide, the pressure rises and first reaches 
the homeostatic pressure pH,h- At this point, tissue H 
stops growing while tissue T continues to proliferate and 
drives the pressure above Phm- As a result, the apop- 
tosis rate of H becomes larger than its division rate, re- 
sulting in its recession. The process continues until tis- 
sue H completely disappears. The winning compartment 
always corresponds to the tissue with the larger homeo- 
static pressure. 

It is interesting to consider the effect of biochemical 
signaling or immunological interactions between the two 



tissues. In particular, consider the case where H resists 
the expansion of T by locally decreasing the homeostatic 
pressure of T. If this decrease is large enough, tissue T 
shrinks and the result of the competition will be reversed 
as compared to the case without signaling. However, if 
we consider a sufficiently large compartment T, the re- 
gion close to H has a negligible contribution to the overall 
compartment growth, and T expands as in the absence 
of signaling. There is a particular size of compartment T 
for which, at the homeostatic pressure pn.hi its average 
growth vanishes: the excess division away from the piston 
exactly balances the excess death close to it. A steady 
state is possible for this particular size of compartment 
T, but it is unstable. This introduces the second impor- 
tant concept of this paper: the existence of a critical size 
beyond which a tumor tends to grow and below which it 
tends to shrink. 

A critical size can also exist due to interfacial ten- 
sion in higher-dimensional geometries, such as the two- 
dimensional organisation of a monolayered epithelium or 
the three-dimensional configuration of a secondary tu- 
mor within the bulk of a host tissue. The concept of tis- 
sue interfacial tension has already been used to explain 
cell sorting of tissues with different adhesive properties 
(Duguay et ai, 2003), and quantified for several tissues 
(Foty, 1996). Tissue interfacial tension can also originate 
from the mechanical contraction of cytoskeletal elements 
at the interface (Lecuit and Lenne, 2007; Schotz et ai, 
2008) . In a spheroid of tissue T located within the bulk of 
tissue H, the excess pressure in T is given by Laplace's 
law: Ap = 27/r, where 7 is the interfacial tension be- 
tween H and T and r is the radius of the spheroid. As a 
result, for small enough radii, the pressure in T is larger 
than pT,h, and T recedes. For large radii however, the 
excess pressure as given by Laplace's law vanishes and 
we recover the previous one-dimensional case where T 
grows. There is again an unstable critical radius Tc for 
which a steady state exists. 

So far, we have considered cell growth and death pro- 
cesses as entirely deterministic, in which case only tumors 
larger than the critical size can grow. However, single 
cells give rise to tumors and metastases (Talmadge and 
Fidler, 1982; Talmadge and Zbar, 1987; Chambers and 
Wilson, 1988). This is possible because cell growth and 
death are stochastic processes. In this paper, we calculate 
the probability for a single cell to give rise to a macro- 
scopic tumor and obtain results that are compatible with 
experimental data on metastatic inefficiency (Luzzi et al, 
1998; Cameron et ai, 2000; Zijlstra et ai, 2002). The 
concepts we use here are similar to those used to de- 
scribe the statistics of nucleation processes as they occur 
in first-order phase transitions. It is well known that nu- 
cleation is easier on surfaces or foreign bodies than in the 
bulk of a system. The same holds true for tumor growth: 
we show that it is more likely for tumors to reach the 
critical size at an interface than in the bulk of a tissue, 
in agreement with experimental and clinical observations 
(Cameron et al., 2000; Weiss, 1985). Hence, in this pa- 
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per, we argue that an unstable critical size for tumor 
growth exists, which is responsible for the inefficiency of 
the metastatic cascade and could account for the pre- 
ferred growth of metastases on surfaces and interfaces. 
We treat only the early stages of tumor and metastatic 
growth, where the heterogeneity of tumors — due to ef- 
fects such as the diffusion of nutrients and growth fac- 
tors or genetic mutations — can be neglected. These ef- 
fects play an important role for larger tumor sizes only 
(Hanahan and Weinberg, 2000). 

Results 

Tissue Rheology and Homeostatsis 

While the notion of homeostatic pressure and density is 
model independent, the details of the tissue dynamics are 
not. Here, we employ a continuous description that we 
expect to be valid for systems large compared to the cell 
size and time scales large compared to the characteristic 
times of individual cellular processes. The local density 
of cells p obeys the continuity equation: 

-p + v ■{pw)^{ka-K)p, (1) 

where v denotes the local velocity of the tissue and 
V-(pv) the divergence of the cell flux. The right hand side 
corresponds to source and sink terms that describe the 
local production and destruction of cells due to cell di- 
vision {kd) and apoptosis (fco). In addition, tissues must 
also satisfy momentum conservation, which, for systems 
where inertia plays a negligible role, reduces to force- 
balance: 

daOap = 0. (2) 

Here, da denotes the partial derivative with respect to 
the coordinate a [a = x,y,z), and summation over re- 
peated indices is implicit; (Tap denotes the total stress 
tensor that we split into a velocity-independent part and 
a dynamic part c'^p- For an isotropic tissue, the velocity- 
independent part reads —pSap, where p is the tissue pres- 
sure discussed above. The viscous part however encodes 
the rheological properties of the tissue in a constitutive 
equation that relates it to the velocity-gradient tensor 
daVfj. Tissues are complex media with a rheological be- 
havior intermediate between those of liquids and solids 
(Foty et al, 1994; Schotz et al, 2008). On timescales 
short compared to their viscoelastic relaxation time, tis- 
sues have a finite shear modulus E of the order of 10^- 
10'' Pa (Forgacs et al, 1998; Engler et al, 2004; Kong 
et al, 2005). For time scales exceeding the largest re- 
laxation time T however, viscoelastic media behave as 
viscous liquids with viscosity rj = Et. Measurements of 
the mechanical response of various cell aggregates sug- 
gest a value of the relaxation time in the range of tens of 
seconds to several minutes (Forgacs et al., 1998; Schotz 



et al. , 2008) , corresponding to a viscosity in the range of 
10^-10^ Pa-s. The fastest division rates of mammalian 
cells are typically of the order of one division per day 
(Weinberg, 2007). Hence, tissue-growth dynamics takes 
place on time scales that are long compared to the char- 
acteristic times of cellular processes, including adhesion 
and detachment of the proteins that insure the integrity 
of the tissue under consideration. Under such conditions, 
it is a general result that the effective rheology on large 
scales appears to be that of a fluid (Frisch et al., 1985). 
We therefore argue that, in the context of tissue-growth 
dynamics, a purely viscous rheology is appropriate, which 
leads to the standard constitutive equation: 

= V {daVf3 + dpVa) ■ (3) 

Under flxed biochemical and biophysical conditions, di- 
vision and apoptosis rates — as well as pressure — are func- 
tions of cell density only. In the absence of a detailed 
knowledge of the pressure and rate dependences as func- 
tions of p, and for the sake of simplicity, we rely on an 
expansion to first order in p — Ph around the homeostatic 
density p^: 

P = X^^iP - Ph)+Ph 
kd-ka = -K{p-ph). (4) 

The parameter x is equivalent to the standard compress- 
ibility of a material and describes the variation of cell 
density with pressure. Similarly, the coefficient k quan- 
tifies how the difference between division and apoptosis 
rates depend on density. Both x ^J^d k are experimen- 
tally accessible parameters that must both be positive to 
insure stability. In Eqs. ([4]), the expansion of the pressure 
p in terms of the cell density p is complementary to the 
expansion of kd — ka known as logistic growth, which is 
a common way to model growth dynamics (Sachs et al, 
2001). As we have stated, the most general dependence 
of the pressure p as well as of the division and apopto- 
sis function kd — ka on the biochemical environment of 
the tissue is encoded in the expansion coefficients x 
K here, which are therefore constants only under fixed 
biochemical conditions. As we focus here on tumors of 
small sizes under steady environmental conditions, such 
a dependence will not be discussed here. But our frame- 
work in principle allows to study more complex situations 
where spatio-temporal inhomogeneities would play a role, 
simply by allowing x ^J^d k to vary. Cell division and 
apoptosis could also be coupled directly to pressure, lead- 
ing to tissue competition as proposed in a similar model 
by B. Shraiman (Shraiman, 2005). Also, studies have 
shown that a defective density sensing of cancerous cells 
can lead to a growth advantage (Chaplain et al., 2006). 
In this paper, we argue that competition for volume is a 
generic property of tissues in mechanical contact, since 
the pressure and effective growth rate functions take the 
form of Eqs. Q close to the steady state density of the 
tissue. 
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Tumor Growth Dynamics 

We are now in a position to show that homeostatic 
regulation intrinsically contains a growth mechanism for 
neoplastic tissues. Consider the growth of a spherical tis- 
sue T located at the center of a spherical compartment 
of finite volume, filled with another tissue H of lower 
homeostatic pressure. The spherical tissue can be either 
a primary tumor developing within the tissue it stems 
from, or a metastasis that has migrated from its original 
location and invaded a foreign organ. The two tissues are 
in mechanical contact, so that the total stress is continu- 
ous at the interface. Eqs. ([iJ-Q must be solved for both 
compartments, taking the location of the interface of the 
two tissues into account. A numerical solution of the as- 
sociated generic growth dynamics is presented in Fig. [2] 
The solution shows that the tissue with higher homeo- 
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FIG. 2: Numerical solution for the cell density (A) and cell 
velocity (B) as functions of space and time during the growth 
of one tissue located in the bulk of another tissue of lower 
homeostatic pressure. Color coding for local cell density and 
velocity is given on the right hand side of panels (A) and 
(B), respectively. Spherical symmetry is assumed. Total inte- 
gration time, compartment size and homeostatic densities of 
tissues T and H are scaled to one. In both plots, the boundary 
between the two tissues is indicated by a black line. Param- 
eters are chosen in order to illustrate the interplay between 
viscous dynamics and compartment growth (see supplemen- 
tary material). 

static pressure grows at the expense of the other one and 
takes over the entire compartment. In vivo however, the 
condition of a fixed finite volume does not hold in general. 
In real tissues, there is often first a displacement of the 
non tumor tissue before anatomical constraints limit the 
total volume available to the system. However, the dev- 
astating effect of malignant tumors stems from the fact 
that they invade and replace the functional tissues. The 
architecture of most tissues leads to a competition for 
volume in the case of neoplastic proliferation. Note that, 
considering the time evolution of the boundary between 
the two tissues only, we get a curve that is reminiscent 
of the well-kown, experimentally-observed Gompertzian 
growth curves (Molski and Konarski, 2003). A quanti- 
tative illustration of this behavior obtained within our 
framework is illustrated in the supplementary material, 
making use of realistic parameters. 

We now examine several effects that can significantly 



alter the tumor growth dynamics as presented above. A 
first example, which is motivated by the structure of 
benign tumors, corresponds to tissue T engulfed in a 
membrane, typically a thin shell of extracellular matrix, 
where the surface tension 7 rises as T expands. If this 
tension increases faster than the radius of T, the addi- 
tional pressure increases and the expansion of T eventu- 
ally stops: a stable steady state exists at r = Tg such that 
PT,h-PH,h = 2-f{rs)/rs, where pT,h and pH.h are the re- 
spective homeostatic pressures of T and H. A numerical 
solution illustrating this case is presented in Fig. |3] This 
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FIG. 3: Numerical solution for the cell density (A) and cell 
velocity (B) as functions of space and time during the growth 
of a tissue engulfed in an elastic membrane, and located in 
the bulk of another tissue of lower homeostatic pressure. As 
in Fig. [2] spherical symmetry is assumed; total integration 
time, compartment size and homeostatic densities are scaled 
to one, and in both plots the boundary between the two tis- 
sues is indicated by a black line. Color coding is similar to the 
one used in Fig. [2] In this solution, the membrane is treated 
as purely elastic and is put under tension above a given ra- 
dius xo = 0.5. The explicit surface tension dependence on the 
boundary location x is given in the supplementary material. 
Additional parameters are given in the supplementary mate- 
rial. (A) The expansion of the inner compartment is similar 
to that of Fig. [2] until, at 2; = 0.5, surface tension begins to 
play a role. The membrane expands until its tension balances 
the pressure difference between the two compartments and a 
stable steady state is reached. (B) Corresponding cell-velocity 
plot. 

dormant state is stable until genetic alterations inducing 
the production of proteases by the tumor cells lead to the 
degradation of the membrane. 

As a second example, we consider the case of a tumor 
that is limited in its growth, for example by nutrient or 
oxygen supply. It is indeed a well-known fact that tumors 
are poorly vascularized before they acquire the capability 
to trigger the growth of new blood vessels via angiogene- 
sis (Folkman and D'Amore, 1996; Weidner, 1991; Hana- 
han and Weinberg, 2000). This limitation has profound 
consequences for their growth dynamics (Preziosi, 2003), 
often leading to the existence of a maximum size of about 
one to two millimeters, where they remain in a "dor- 
mant state" until the induction of angiogenesis (Folk- 
man and D'Amore, 1996; Weinberg, 2007). In Fig. |4] 
we present a numerical integration of the growth dynam- 
ics of a nutrient-limited tumor in the bulk of a healthy, 
well vascularized tissue. While we assume a homoge- 
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FIG. 4: Numerical solution of the growth dynamics with 
the geometrical arrangement of Fig. [2] when growth rates are 
nutrient-limited. Nutrients diffuse into the inner compart- 
ment through the tissue interface (see supplementary mate- 
rial). (A) The inner compartment starts growing as in the 
case of Fig. [2] but asymptotically reaches a maximum size. 
(B) Cells proliferate at the surface of the inner compartment 
which is rich in nutrients and die at the center where nutrients 
are scarce, resulting in an inward flow of cells. Color coding 
is similar to the one used in Fig. |2] and [S] but it now allows 
for negative values required by the inward flow. 

neous and high enough concentration of nutrients for the 
healthy tissue, the neoplastic tissue is only supplied with 
nutrients via diffusion through its surface. Since nutri- 
ent diffusion is fast compared to growth dynamics, we 
calculate the nutrient concentration profile by solving a 
steady-state diffusion equation, taking into account nu- 
trient consumption due to cell metabolism and division 
(see supplementary material). We choose functional de- 
pendences of the division and apoptosis rates on the nu- 
trient concentration and cell density that correspond to a 
biological behavior: under very low concentrations of nu- 
trients or oxygen, cells tend to die, but a limited supply of 
nutrients can also decrease cell division by triggering cell 
differentiation, inducing a quiescent cell state or favoring 
adaptation of the metabolism of the cells to the new en- 
vironment. In agreement with what is known about the 
internal structure of dormant tumors, cells divide at the 
boundary where they get enough nutrients, and die at 
the center. This creates a steady state flow of cells from 
the surface toward the center of the tumor and thereby 
a constant cell turnover that is favorable to mutations. 

Note that the typical length-scale at which the diffu- 
sion of nutrients becomes a limiting factor is of the order 
of millimeters, the size of a dormant tumor (Folkman and 
D'Amore, 1996; Weinberg, 2007). This scale is very large 
compared to the size we estimate for the critical radius 
introduced above. In the following, when considering the 
nucleation process of micro-tumors, we therefore assume 
a homogeneous, high enough concentration of nutrients. 

Critical Size and Stochastic Growth Dynamics 

A third effect that can modify tumor growth dynamics 
is the presence of a constant tissue interfacial tension as 
introduced above. When H and T are at their home- 
ostatic densities, there exists a particular radius Tc at 



which mechanical equilibrium is reached, but this equilib- 
rium, given by px.h — PH,h = 27/rc, is unstable. Numer- 
ical solutions illustrating the growth dynamics around 
this critical radius are presented in Fig. |5] 




Time 

FIG. 5: Tissue boundary as a function of time during the 
growth of a tissue located in the bulk of another tissue of 
lower homeostatic pressure, with interfacial tension and in 
spherical geometry. Parameters are given in the supplemen- 
tary material, together with Vc = 0.5 and 7=1. The different 
curves show the dynamics for the following initial values ro of 
r: 0.49, 0.499, 0.4999, 0.501 and 0.51. Both tissues start out 
at their homeostatic densities. 

Given the existence of such an unstable critical radius, 
the question arises as to how a metastasis — or a primary 
tumor — can grow within a healthy tissue since in gen- 
eral it originates from a single cell (Talmadge and Fidler, 
1982; Talmadge and Zbar, 1987; Chambers and Wilson, 
1988). The answer stems from stochasticity, an aspect 
of the dynamics that has been ignored in the descrip- 
tion so far. The importance of stochasticity in growth 
processes has already been recognized in various situa- 
tions (Nowak et ai, 2003). Under the assumption that 
stochastic tumor growth is a Poisson process, the evo- 
lution of the probability P{n, t) for a spherical tumor T 
inside a healthy tissue H to contain n cells at time t can 
be described by a master equation: 

+r-+iP(n+l,i), (5) 

where = nkd and = nka are the rates at which a 
tumor grows or shrinks from n to (n -|- 1) or (n — 1) cells, 
respectively. The rates kd and ka depend on n, and we 
model their dependence in the following way: For tumors 
small compared to the size of the healthy compartment, 
the healthy tissue is only slightly perturbed away from its 
homeostatic state. Thus, the pressure inside the tumor 
is given by Laplace's law: = PH,h + 27/r. Therefore, 
the division and apoptosis rates of a spherical tumor of 
radius r are given by: 

kd/a = -Kd/a XT (pH,h + — - PTJij + ^0- (6) 
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Here, Kd, Ha and fco are three phenomenological coef- 
ficients that enter the linear expansions of kd and ka, 
similarly to k in Eq. ^. To ensure the proper behavior 
as a function of the cell density p, Hd needs to be posi- 
tive and Ka negative. Both equations for kd and ka share 
the same constant fco such that Eq. Q is satisfied with 

In the master equation ([5|, the rates 
and r~ are then given to leading order by: 



-Krf/a XT PH,h + 27 



3?^ 



1/3 



PT,h 



(7) 



and fco > fixes the amount of cell turnover — and 
thereby the amount of stochasticity — in the system. 

For an analytic treatment, we map this growth process 
onto a random walk with sinks at n = and n = nmax, 
which results in a linear birth-death process where all 
tumors either disappear or reach macroscopic sizes when 
time goes to infinity. The so-called "splitting probability" 
lig — namely the probability for a single cell to reach the 
size n — nmax and not disappear in the lower sink n — — 
is given by (Van Kampen, 2007): 



^9 = 



1 



-1 r. 



(8) 



In Fig. |6] we present this analytic result together with 
the results of a Monte Carlo simulation of Eqs. ([5])-([7]) 
based on a Gillespie algorithm (Gillespie, 1977). 

Finally, it is interesting to consider the growth of neo- 
plastic semi-spheroids on tissue boundaries. The internal 
pressure of a semi-spheroid is again given by Laplace's 
law, but for the same radius, the number of cells within 
the tumor is half that of a complete spheroid. Therefore, 
the critical number of cells of a tumor at an interface is 
only half that of a tumor in the bulk. We can directly 
read the resulting effect on the growth probability out of 
Fig.|6] For example, between the critical sizes of five and 
ten cells, we obtain a growth probability ratio of 8.6. This 
characterizes a significant preference for tumors to grow 
on surfaces, an effect that has been observed experimen- 
tally (Cameron et al, 2000) and clinically (Weiss, 1985). 
Note that other mechanisms — such as the adhesion of 
cancerous cells to the extracellular matrix, which leads to 
a different initial distribution of metastatic cells — might 
also play a role. 



Discussion 

In this work, we have shown that homeostatic regu- 
lation of cell density and pressure leads to a competi- 
tion for space between tissues in mechanical contact. We 
have proposed that an increased homeostatic pressure is 
a characteristic trait of tumors. Gedanken experiments 
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FIG. 6: Splitting probability as a function of the critical cell 
number. The black dots are the fraction of tumors with a 
size larger than the critical radius after 14 days of evolution, 
starting from a single cell in a Monte Carlo simulation (see 
Eqs. |[5| and ([7|). The black line represents the analytic result 
for the splitting probability as given by Eq. Q. Parameters 
are chosen such that far above the critical radius the division 
rate reaches a value of one division per day and the apoptosis 
rate is negligible, fco in Eq. ([7| is chosen to be 0.9 division per 
day. Note that the probability for growth on the compartment 
surface can be read out from the plot using half of the critical 
cell number of the process in the bulk. Upper inset: Fraction 
of tumors that are above the critical size. This fraction relaxes 
to a constant value in a Monte Carlo simulation. Lower inset: 
Distribution of tumor sizes in a Monte Carlo simulation after 
10 days of evolution of 10^ individual tumor cells. 



for measuring the homeostatic pressure and growth rates 
have been discussed and could be realized using current 
experimental techniques. For example, Helminger et al. 
have shown that tumors growing in agarose gels prolif- 
erate until the pressure exerted by the gel reaches 45- 
120 mmHg (Helminger et al, 1997). This gives an esti- 
mate of the homeostatic pressure defined in this paper. 
Such numbers are compatible with the pressure typically 
generated by actin polymerization (Footer et al, 2007; 
Marcy et al., 2004). An experiment of particular interest 
that has not been performed so far is the competition 
between a healthy and a tumor-like tissue separated by a 
piston, which could prove that mechanical effects are im- 
portant for tumor growth. If the growth mechanism dis- 
cussed in this paper is relevant to tumor growth, it would 
be interesting to measure the homeostatic pressures char- 
acterizing healthy and neoplastic tissues, together with 
their dependences on their biochemical environments. Of 
particular interest would be their dependence on oxygen, 
nutrients, growth factors and drugs. 

The second concept introduced in this paper is the exis- 
tence of a critical size for tumor growth due to biochemi- 
cal, immunological or mechanical surface effects that can 
outbalance the bulk growth advantage of the neoplas- 
tic tissue for small tumor sizes. We show that this in- 
teraction can be responsible for the inefficiency of the 
metastatic cascade after extravasation. The growth of 
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very few metastases to macroscopic sizes cannot be ex- 
plained by size-independent growth rates, which yield a 
probability distribution of metastatic cell clusters that 
decays exponentially with cluster size. Instead, the exis- 
tence of a critical size yields realistic values for metastatic 
inefficiency and a distribution of tumor sizes compati- 
ble with experimental observations (Luzzi et al., 1998; 
Cameron et al, 2000; Zijlstra et ai, 2002) (see Fig. ^. 
Fig. [6] also shows the dependence of metastatic ineffi- 
ciency on the critical size: a small change in tissue- 
tumor interaction such as an increased interfacial ten- 
sion can dramatically lower the probability for macro- 
scopic growth. As an illustration of this effect, consider 
a metastatic tissue with a critical cell number of 5 in a 
given environment. Let us compare this situation with 
that of the same tissue placed in another environment 
where its interfacial tension is now twice as large, a situ- 
ation that is well within natural variations (Foty, 1996). 
While for the first environment, about 3 in 100 metastatic 
cells form a macroscopic tumor, in the second environ- 
ment, with a critical cell number of 40, less than 2 in ten 
million manage to do so. This corresponds to a differ- 
ence in metastatic efficiency of five orders of magnitude. 
We propose that this effect could account for the strong 
tissue specificity of metastatic growth that underlies the 
"seed and soil hypothesis" (Weinberg, 2007). 



The concept of homeostatic pressure presented here 
is not an alternative to the cellular and genetic mecha- 
nisms involved in tumor growth, but rather a different 
level of description. Indeed, we propose that some of the 
fundamental biological deregulations that are character- 
istic of neoplastic cells lead to an increased homeostatic 
pressure. The framework presented here can be used to 
explicitly take into account such well-known properties. 
It can also be generalized to incorporate more general 
features of biological tissue behavior. For example, on 
long time scales, genetic instability as well as senescence 
render tissue properties time dependent. This could be 
incorporated into our framework using techniques sim- 
ilar to those of Hallatschek et al. (Hallatschek et al, 
2007; Hallatschek and Nelson, 2008), as well as those of 
multiscale models of tumor growth (Ribba et al., 2006; 
Macklin and Lowengrub, 2007; Wise et al, 2008). 
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Supplemental Material 



I. SPHERICAL GROWTH 

As an example of a tissue growth competition for which 
we can solve the complete dynamics given by the Eqs. 1- 
4 of the main text, we examine the growth of a spherical 
tissue located in the center of a spherical compartment 
filled with another tissue of lower homeostatic pressure 
and enclosed by a rigid boundary. For now, we neglect 
all surface tension effects. The force balance condition 
(Eq. 2, main text) in three dimensions takes the form: 



dr 



2{a'„ - (j'gg) 



dp 
dr 



(9) 



where r, 9 and are the spherical coordinates, while the 
constitutive equation (Eq. 3, main text) reads: 



2r? 
2r? 



dv 
dr 

V 



(10) 



The continuity equation (Eq. 1, main text), together with 
the expansion of kd — ka to first order in p — ph (Eq. 4, 
main text), gives: 



dp 
di 



1 d 



(11) 



These equations need to be solved for the whole system 
composed of the two tissues, together with the mov- 
ing boundary between them. Boundary conditions are 
composed of two parts: (a) in the center and at the 
rigid external wall, the velocity field vanishes, such that 
w(r = 0) = w(r = i?) = 0; (b) at the interface of the two 
tissues, the velocity field and the stress tensor are contin- 
uous. The continuity of the velocity field at the interface 
leads to the following equation for the time-dependent 
location x{t) of the tissue boundary: 



dx 
dt 



v(x). 



(12) 



The force-balance condition ([9| can be integrated us- 



ing: 



dr 



-{a'„ 



from the constitutive equation (10 1 to give 



2r, 



dv 2v 
dr r 



(13) 



(14) 



In the absence of surface tension, the integration con- 
stant pcxt is the external pressure imposed by the rigid 



wall to satisfy the boundary condition of vanishing ve- 
locity. The growth dynamics consisting of Eqs. (11) and 
( 14 1 for each of the two compartments — together with the 



moving boundary condition Eq. ( 12 1 — can be solved nu- 
merically using a finite-difference method (Press et al, 
1992). Results with the parameters of Table |l] are dis- 
played in Fig. 2 of the main text and show how the inner 
tissue takes over the whole compartment. 

In the main text, a constant interfacial tension is in- 
troduced between the two tissues. It is shown that this 
effect leads to the existence of an unstable critical ra- 
dius in spherical geometry. The growth dynamics with 
an unstable critical radius is illustrated in Fig. 5 of the 
main text. However, biologically relevant situations may 
involve tissues enclosed in membranes whose tensions in- 
crease as the inner tissue grows. This is for example the 
case for some benign tumors that undergo growth ar- 
rest due to the extracellular membrane engulfing them. 
In that case, the surface tension is now dependent on the 
location x of the boundary between the two tissues. For a 
purely elastic membrane that is put under tension above 
a given radius xq, we have: 

l{x) = jo^-^^^0{^-xo), (15) 

Xq 

where 0(x) is the heaviside step function. A numerical 
solution of the growth dynamics with this type of surface 
tension is presented in Fig. 3 of the main text (parame- 
ters are given in Table |T] together with 70 = 5 in scaled 
units). 



II. TUMOR GROWTH DYNAMICS WITH 
REALISTIC PARAMETERS 

The growth rate of tumor cells is possibly very slow 
compared to the viscous relaxation time. In such a case, 
it can be assumed that the cell density in each compart- 
ment is constant. In the absence of surface tension, the 
pressures in the two compartments balance, leading to 
an equation relating the two densities: 

Xh^{PH ~ PH.h)+PHM = Xt^{PT - PT,h)+PT,h- (16) 

The change in density in each compartment has contri- 
butions coming from the total cell division and apoptosis 
taking place in the compartment, as well as from the 
movement of the boundary x(t): 



Pt 
P'h 



-kt {pt - PTji) Pt 



i- Pt 

X 



X 



-kh {ph - PH,h) Ph + 3 ^ 3 



X 

X R?' — X'' 



PH. (17) 



This system of differential equations can be solved nu- 
merically. Results with the parameters given in Table [IT] 
are displayed in Fig. [7] 
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FIG. 7: Radius and cell-density as functions of time for a 
growing tumor in conditions identical to those of Fig. 2 of the 
main text, but with estimates of realistic parameters. Cell 
density is one per 10^ /im^, and cell division rate is one per 
day. The homeostatic pressure difference and the tissue vis- 
cosity are 10'^ Pa and 10* Pa-s, respectively (Helminger et ai, 
1997; Kruse et al, 2000; Forgacs et al, 1998). Tissue com- 
pressibility is 10~^ Pa~^-m~^ (Tschumperlin et al, 2004). 
With these parameters, a separation of timescales between 
the slow cell division and the comparatively fast viscous dy- 
namics occurs, such that we can assume homogeneous cell 
densities in both compartments. Starting from a single cell 
and neglecting nutrient coupling, a tumor needs about 100 
days to fill a compartment with radius lOOpm. 



III. NUTRIENT-LIMITED GROWTH 

When studying the nutrient-limited growth of a tu- 
mor prior to angiogenesis, the dependences of the cell 
division and apoptosis rates kd and ka on the nutrient 
and cell densities p„ and p are constructed using two 
assumptions: (a) below a given concentration of nutri- 
ents per cell ci, cells stop dividing; (b) below a second, 
lower concentration of nutrients per cell C2 , cells undergo 
apoptosis. We model these properties with the following 
functions |3] : 




A 

0.75 
0.5 
0.25 



/On 



FIG. 8: The cell division rate kd (A) and the negative apop- 
tosis rate —ka (B) as functions of the cell density p and the 
nutrient concentration p„, as given by Eq. (181. Parameters 
are given in Table IITI] 



we illustrate the dependence of this combination on cell 
density and nutrient concentration in Figs. [9]and [T0| 
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FIG. 9: Difference kd — ka as a function of cell density and 



nutrient concentration, as given by Eq. (18 1. Parameters are 
given in Table |III[ (A) Top view in grey-scale coding. (B) 
Three-dimensional view. 



1 + exp (a(p -ph + Ap)) I (^_/3^(&L _ ci)) 

ka — 

ki I h 

1 + exp {-a{p -ph-Ap)) C2] 

(18) 

Here, ki tunes the amplitude of cell division and apopto- 
sis in the system as functions of cell density, as ^2 tunes 
how strongly cells die when deprived of nutrients. The 
parameter a tunes how sharply cell start to die or prolif- 
erate as the homeostatic density is passed. It is the same 
in both functions kd and ka, such that fc^ — = at 
P — Ph for large concentrations of nutrients. Ap sets the 
amount of cell turnover at homeostatic density. Finally, 
/3i and (32 tune how sharply the cell division and apopto- 
sis rates change as the critical concentrations of nutrients 
per cell ci and C2 are passed. We illustrate the depen- 
dence of kd and ka on p and p„ in Fig. |8] with param- 
eters given in Tables |T] and Since only the difference 
kd — ka enters the growth dynamics (Eq. 1-4, main text). 
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FIG. 10: Effect of nutrient concentration on the difference 
kd — ka as. & function of cell density. Parameters are given in 
Table [nil ^t high nutrient concentrations, regulation towards 
the homeostatic density = 1 is intact, while at low nutri- 
ent concentrations cell division drops significantly. At very 
low nutrient concentrations, cells undergo apoptosis at a high 
rate. 

To describe how nutrients are distributed in the sys- 
tem, we suppose that they diffuse freely throughout the 
system with a given diffusion constant while being 
consumed by living cells for their metabolism and their 
growth. Metabolism uptake happens at a given rate ppn 
that is proportional to the available concentration of nu- 
trients, and growth dependence is described via an extra 
consumption term proportional to the number of cell di- 
vision kdp with a coupling constant A. For simplicity, we 
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suppose no effect of cell apoptosis on nutrient uptake. We 
finally suppose that nutrient diffusion is very fast com- 
pared to tissue growth, such that only the steady-state 
diffusion equation needs to be considered: 



in addition to the lower sink at n„ 



0, all clus- 



1 d_ 
dr 



jdpn 

dr 



XkdP - fipnP = 0. 



(19) 



Boundary conditions are as follows: the nutrient con- 
centration is homogeneous and constant in the healthy 
compartment and the flow of nutrients vanishes at the 
center of the tumor compartment. 

To compute the growth dynamics of a spherical tumor 
coupled to nutrient diffusion through its surface, we use 
the same method as in Section I, while solving the steady- 
state diffusion equation ( 19 ) at every timestep. The re- 



sult with the parameters given in Table III is shown in 
Fig. 4 of the main text. The characteristic nutrient pro- 
file in a tumor that is nutrient-limited in its growth is 
given in Fig. |11[ 




FIG. 11: Typical nutrient profile as a function of time in 
a growing tumor that arrests in a dormant state. The flux 
at the center vanishes and the outside compartment has a 
constant concentration. Parameters correspond to those given 
in Tables m and Iml 



IV. STOCHASTIC DYNAMICS 

We solve the master equation (Eq. 5, main text) to- 
gether with the rates given by Eq. 7 (main text), both 
analytically and numerically. Imposing an upper sink at 



ters of cells end up in one of the two sinks when time 
goes to infinity. The analytic solution — given by Eq. 8 
(main text) — gives the splitting probability, namely the 
probability for a cluster originating from a single cell to 
reach the upper sink 7i,nax when time goes to infinity. Nu- 
merically, we use a Monte Carlo simulation of the master 
equation 5 (main text) based on the Gillespie algorithm 
(Gillespie, 1976), which gives information on the tempo- 
ral evolution of the process leading to the analytic result 
Eq. 8 (main text) for long evolution times. 

The Monte Carlo simulation is implemented accord- 
ing to the following standard procedure: at each step 
of the growth process (corresponding to a new division 
or apoptosis event), the algorithm first sets the time de- 
lay between this event and the previous one, and then 
chooses whether apoptosis or division takes place. These 
two choices are made by generating two random numbers 
xi and X2 in the interval [0, 1] using a Mersenne Twister 
algorithm (Matsumoto and Nishimura, 1998). Xi deter- 
mines the time delay St between two events as: 



6t=- 



log(a::i 



(20) 



where and r~ are the growth and death rates given by 
Eq. 7 (main text). The stochastic variable X2 determines 
whether growth or recession takes place: growth is chosen 



if [(r 



')x2] < r^, recession otherwise. 



Parameters are chosen as follows: the upper sink is at 
f^max = 10^ cells. The parameters that enter the expres- 
sion for the rates (Eq. 7, main text) are such that, at very 
large radii, the tumor divides at a rate of one division per 
day on average, while having a vanishing probability to 
shrink. This yields: 



l^d/a XT 



k 



+/- 



with 



PT,h - PH,h 



lim -^ = 1 

n— »oo n 



(21) 



hm 

n—^OQ Jl 







(22) 



in invert units of days. Here, ko is an adjustable parame- 
ter in the interval [k~^^, k^_^^] that tunes the amount 
of stochasticity in the system. Finally, we impose a crit- 
ical radius Tc corresponding to a critical number of cells 
Tie at density PT,h, and choose the surface tension accord- 
ingly as: 



PTM - PH,h ( 371c 



\A.TipTM 



1/3 



(23) 



The parameters used in Fig. 6 of the main text are given 
in Table HVI 
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Tables 



PaiaiiK?t(?r 


\a.luc 


Dcscriptitju 


Xo 


0.01 


initial interface location 


Po 


1 


initial density 


vo 





homogeneous initial velocity 


ievol. 


650 


total time of evolution in Fig. 2 (main text) 




600 


total time of evolution in Figs. 3 and 5 (main text) 




3000 


total time of evolution in Fig. 4 (main text) and Fig. 5 here 


K 


10 


division constant (both tissues) 


X 


0.2 


compressibility (both tissues) 


V 


50 


viscosity (both tissues) 



TABLE I: Parameters used to compute the growth dynamics of spherical tumors in Figs. 2, 3, 4 and 5 of the main text and 
Fig. 5 here. Parameters are scaled in units of the total compartment size R, the homeostatic pressure difference Ap between 
the two tissues and the homeostatic densities ph which we assume to be identical for the two tissues. Total times are scaled to 
one in the figures. 



Parameter 


Value 


Description 


R 

Xo 

Ph 
Ap 

X 

K-x- 
V 


100 urn 
10 /um 
0.001 A«m^ 
1000 Pa 

1 d-i 
10"* Pa-s 


total compartment radius (Weinberg, 2007) 
initial interface location 

homeostatic density (both tissues) (Weinberg, 2007) 

difference in homeostatic pressures (both tissues) (Helminger et al., 1997) 
compressibility (both tissues) (Tschumperlin et al., 2004) 
maximum division rate (both tissues) (Weinberg, 2007) 
tissue viscosity (Forgacs et al., 1998) 



TABLE II: List of parameters used to compute the growth dynamics of spherical tumors in Fig. 1. 
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Parameter 


Value 


Description 


Cl 





nutrients per cell for induction of apoptosis 


C2 


0.6 


nutrients per cell for arrest of proliferation 




10 


maximum cell division and apoptosis at high nutrient concentration 


h 


500 


apoptosis rate coefficient at starvation 


a 


10 


response coefficient of cell division and apoptosis to cell density 


I3i 


50 


response coefficient of cell division to nutrient concentration per cell 




500 


response coefficient of cell apoptosis to nutrient concentration per cell 


Ap 


0.1 


shift in cell division and apoptosis tuning the amount of 
cell turnover at homeostatic density 


D„ 


1 


nutrient diffusion constant 


A 


6 


nutrient consumption for proliferation 


M 


1000 


nutrient consumption for metabolism 



TABLE 111: Parameters used to generate the plots of Fig. 4 of the main text and Figs. 2-5 here in the coupling of the growth 
dynamics to nutrients as they enter in Eq. ( |18[ ). 



Parameter 


Value 


Description 


k+ 

"'n — >-oo 
k~ 

""n — »oo 


1 d-i 
d-i 
0.9 d-^ 


n infinity value of /n 
n infinity value of /n 
cell turnover at homeostatic density 



TABLE IV: Parameters used in Fig. 6 of the main text. 



